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Abstract 

We compute charm and bottom quark masses in the quenched approximation and 
in the continuum limit of lattice QCD. We make use of a step scaling method, previ- 
ously introduced to deal with two scale problems, that allows to take the continuum 
limit of the lattice data. We determine the RGI quark masses and make the con- 
nection to the MS scheme. The continuum extrapolation gives us a value m^ GI = 
6.73(16) GeV for the 6-quark and rn^ 1 = 1.681(36) GeV for the c-quark, corre- 
sponding respectively to mf^m/P) = 4.33(10) GeV and mf^(mf^) = 1.319(28) 
GeV. The latter result, in agreement with current estimates, is for us a check of the 
method. Using our results on the heavy quark masses we compute the mass of the 
B c meson, M Bc = 6.46(15) GeV. 
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1 Introduction 



Quark masses are fundamental parameters of the QCD Lagrangian. Their 
accurate knowledge is required in order to give quantitative predictions of 
fundamental processes. A direct experimental measurement of quark masses 
is not possible because of confinement, and their determination can only be 
inferred from a theoretical understanding of the hadron phenomenology. The 
calculations can be numerically performed through different strategies, de- 
pending upon the quark flavor and the available computational facilities. Ac- 
curate non-perturbative measurements of the u, d, s and c quark masses have 
been obtained from a straightforward comparison of hadron spectroscopy and 
lattice QCD predictions (see [1] for a recent review). 



The situation is different for the b quark, with present computers capa- 
bilities. In principle, the b quark mass could be extracted, similarly to lighter 
flavors, by looking at the heavy-light and heavy-heavy meson spectrum. How- 
ever, heavy-light mesons are characterized by the presence of two different 
scales, i.e. Aq C d, that sets the wavelengths of the light quark, and the heavy 
6-quark mass. Managing these two scales in a naive way would require a very 
large lattice. Indeed, this should contain enough points to properly resolve the 
propagation of the heavy quark and to make the light quark insensitive to 
finite volumes effects. A typical size would be O(100 4 ) points, hardly afford- 
able in terms of current memory and CPU time. The heavy-heavy mesons 
case is simpler, being characterized by a single scale. However, the heaviness 
of the bound state makes the exponential decay of the meson correlation func- 
tions too fast and the ground state effective mass, at large time separations, 
cannot be disentangled from the numerical noise (at least on single precision 
architectures) . 

The task of determining the 6-quark mass has been faced in literature by 
resorting to some approximations of the full theory, e.g. HQET on the lattice 
[2], lattice NRQCD [3], or QCD sum rules [4,5]. A novel approach recently 
introduced, based on the non-perturbative renormalization of the static theory 
and its matching to QCD, has lead to a very precise determination of the b- 
mass in the static approximation [6,7]. 

An alternative approach to the bottom quark physics, based on finite 
size scaling, has been proposed in a previous paper [8], where it has been 
applied to the heavy-light meson decay constants. The main advantages of 
this step scaling method (SSM) are that the entire computation is performed 
with the relativistic QCD Lagrangian and that the continuum limit can be 
taken, avoiding the unfeasible direct calculation. In order to implement the 
SSM, a finite size scheme is required, and we adopt the Schrodinger Functional 
(SF) as the most useful framework. This paper is devoted to apply the SSM 
to the study of the heavy meson spectrum in order to obtain from it the first 
determination of the 6-mass in the continuum limit of lattice regularization 
and quenched approximation. 

The paper is organized as follows: section 2 introduces the main ideas 
underlying the SSM; in section 3 details are given on its specific implementa- 
tion through the SF. In section 4 we discuss the predictions of HQET on the 
heavy-light step scaling functions. Section 5 includes the analysis and numer- 
ical results for each step. Section 6 contains the results on a physical volume 
and in the continuum limit. The conclusions are drawn in section 7. 
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2 The step scaling method 



The SSM has been designed in order to deal with two scale problems in 
lattice QCD [8], and it has been shown to work successfully in a first estimate 
of the _B-meson decay constant at finite lattice spacing. A detailed explanation 
of the method can also be found in [9] , where the calculation of the _B-meson 
decay constant has been performed in the continuum limit. Here the SSM is 
reviewed to set the notation and explain why it can be used also in the case 
of the 6-quark mass. 

On very general grounds, the method can be conveniently used in order 
to compute physical observables 0(E e , E h ) depending upon two largely sepa- 
rated energy scales E e and E h (E t <C E h ), where direct simulation, without 
introducing big lattice artifacts, would require a very demanding computa- 
tional effort. The main assumption of the method is that the finite size effects 
affecting O have a mild dependence upon variations of the high energy scale 
and are controllable from a numerical point of view. Finite size effects can be 
obtained from the ratio, o"o, of the observable O computed on two different 
finite volumes, e.g. L and 2L, 

°° (E " Eh ' L)= 0(E t ,E h )\ L (1) 

The step scaling function oq should simplify in the region where Eh Eg. In 
principle, a total decoupling of Eh would determine an absolute insensitivity 
to variations of this scale 

a (E e , E h , L) ~ a (E e , L), E h ^> E e (2) 

In practice, Eh never completely decouples, but the mild residual dependence 
can be suitably parametrized. A typical situation is when the residual depen- 
dence upon E h is linear in 1/E h [8]. In the case of a flavored meson where E h 
is the heavy quark mass, this amounts to state that a heavy quark effective 
theory expansion is valid on finite volumes. Hence, eq. (2) has to be corrected 
as 

rr (F F T) - a (F M I , ^\^L) , ^ 
&o{Ee, E h , L) - a {E e , L) -\ 1 h . . . {A) 

Eh E h 

The number of terms that have to be taken into account in the numerical 
calculations depend upon the particular observable (see sec. 5.2). The ansatz 
(3), for the particular case of the heavy-light meson masses it is supported 
by HQET, it has been numerically checked and allows to extrapolate the 
knowledge of the step scaling function to regions of phenomenological interest 
where direct simulations are too expensive. The computation of the observable 
O proceeds according to the following lines. 
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First, O is computed on a small finite volume Lq, where the high energy 
scale Eh can match its phenomenological value with the lattice cutoff much 
larger than E h . This computation is clearly unphysical, because the fmiteness 
of the volume produces a distortion of the result, that cannot be compared 
with the experimental value. 

Second, the step scaling function is used in order to evolve this finite size 
measurement to a larger volume, according to 

0(E e , E h , Loo) = 0(E e , E h , L Q ) a (E e , E h , L ) a D (E e , E h , 2L ) ... (4) 

In the large volumes the step scaling functions are evaluated at the high energy 
scale Eh by extrapolation, relying on the parametrization of eq. (3). 

Each step of the calculation can be extrapolated to the continuum limit. In 
order to match subsequent steps, a non-perturbative knowledge of the lattice 
spacing a(g ) as function of the bare coupling is required [10,11]; the range of 
validity of the non-perturbative parametrization of a(g ) has been extended to 
very small couplings by a renormalization group analysis in [12]. Throughout 
the paper we use r = 0.5 fm. 

The above strategy has been straightforwardly applied in this work to the 
case of the heavy-light meson masses, where the observable O is the meson 
mass M, extracted from non-perturbative simulations of two point correlators, 
and the energy scales E h and Eg are the heavy and the light quark masses 
[mh^rni). The light quark mass is extrapolated to values corresponding to 
the physical u and s masses, while the heavy quark mass is determined by 
comparing the meson mass to the experimental spectrum. The heavy-heavy 
mesons case should not present major finite size effects, being characterized 
by a single heavy scale. However, as we will see, the application of the SSM 
method allows to solve the problem of the fast decay of the meson correlation 
functions. 

3 Renormalization 

The step scaling function is calculated within the SF [13,14], which has al- 
ready been applied to a number of different finite size problems [15,16,17,18,7]. 
In particular, we use a topology T x L 3 with periodic boundary conditions on 
the space directions, Dirichlet boundary conditions along time and the follow- 
ing set of parameters 

T = 2L, C = C' = 0, 6 = (5) 

Here C and C represent the boundary gauge fields and 8 is a topological phase 
which affects the periodicity of the fermion boundary conditions. Lattice dis- 
cretization is performed using non-perturbative 0(a) improved clover action 
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[19] and operators. In order to set the notation, let 



A ll (x)=il> i (x)'y li 'y 5 ip j (x) 
P(x)=ip i (x)^ j (x) 
V li (x)=i/> i (x)'y IM il> j (x) 

T fiu (x)=ip i {x)- f ^ u ilj j {x) (6) 

be the axial current, the axial density, the local vector current and the tensor 
bilinear operator respectively (i and j are flavor indices). The improvement of 
the axial and vector currents is obtained through the relations 

A I li (x)=A IM (x)+ac A d^P(x) 

V^(x) = V fl (x) + ac v d v T^(x) (7) 

where <9 M = (d^ + d*)/2 and <9 M , d* are the usual forward and backward lattice 
derivatives respectively. For what concerns the improvement coefficients, in 
the case of c A we use the non-perturbative results of [19] while, in the case of 
cy, we use the non-perturbative data only at the values of the bare couplings 
where they exist [20] and perturbative results otherwise [21]. The correlation 
functions used to compute pseudoscalar and vector meson masses are defined 
by probing the previous operators with appropriate boundary quark sources 
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/pW = 4E<C,(y)^W> (8) 



2 y,z 

/^o) = -^E<Ci(y)7*Ci(«)^)> 

y,z 

where Ci(y) an d Ci(y) can De considered as quark and anti-quark boundary 
states. 

The so-called bare current quark masses are defined through the lattice 
version of the PCAC relation 



wi _ go/ApEo) + ac A dld Q fp(x ) 
2f P (x ) 



These masses are connected to the renormalization group invariant (RGI) 
quark masses, according to the definitions given in [22], through a renormal- 
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ization factor which has been computed non-perturbatively in [16]: 

ami + arrij ' 



mf 1 = Z M (g ) 



l + (b A - bp) 



< 7 (So) 



(10) 



where am,- is defined as 



ami 



1 



The combination b A — bp of the improvement coefficients of the axial current 
and pseudoscalar density has been non-perturbatively computed in [23,24]. 
The factor Zm(qo) is known with very high precision in a range of inverse 
bare couplings that does not cover all the values of f3 used in our simulations. 
We have used the results reported in table (6) of ref. [16] to parametrize 
Zm{9q) in the enlarged range of f3 values (5.9, 7.6). 

The RGI mass of a given quark is obtained from eq. (10) using the diagonal 
correlations 



mf GI 



m? GI 



(12) 



From non-diagonal correlations in eq. (10) one obtains different 0(a) improved 
definitions of the RGI i-quark mass for different choices of the j-flavor: 

m™ = 2m™ 1 - mf GI (13) 



All these definitions must have the same continuum limit because the depen- 
dence upon the j-flavor is only a lattice artifact. Further, for each definition 
we use in eq. (9) either standard lattice time derivatives as well as improved 
ones [23,24]. 

Another non-perturbative 0(a) improved definition of the RGI quark 
masses can be obtained starting from the bare quark mass 

mf GI = Z M (g ) Z(g ) [1 + b m am*] m, (14) 
where the improvement coefficient b m and the renormalization constant 

Z(9o) = ^ (15) 



have been non-perturbatively computed in ref. [23,24]. 

Equations (12), (13) and (14) give us different possibilities to identify the 
valence quarks inside a given meson (fixed by the values of the bare quark 
masses). The procedure is well defined on small volumes because the RGI 
quark mass is a physical quantity that does not depend upon the scale, given 
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Fig. 1. The figure shows the excited states finite size effect as predicted by a simplified 
qualitative model discussed in the text. The effective mass step scaling functions are 
plotted, at different values of the meson masses, as functions of the volume. 

in the SF scheme by the volume, and is defined in terms of local correlations 
that do not suffer finite volume effects. Each pair (mf GI ,m^ GI ^j fixed a priori 
is matched, changing the values of the hopping parameters, by the different 
definitions of equations (12), (13) and (14), and leads to values of the corre- 
sponding meson mass differing by 0(a 2 ) lattice artifacts. We take advantage 
of this plethora of definitions by constraining in a single fit the continuum 
extrapolations (see Fig. [2,4,7]). 

The meson masses are extracted from the so-called effective mass 



where fx is one of the correlations defined in (8). On a physical volume, 
this quantity exhibits a plateau in the time region where the ground state 
dominates the correlation and no boundary effects are present. On a small 
volume the effective mass is affected by two different finite volume effects: the 
first one is due to the compression of the low energy wavelength (if present) 
and the second one comes from the presence of the excited states contribution 
to the correlation (no plateau). A good definition of the meson mass that suits 
the step scaling method is obtained by choosing the value of the effective mass 
at xq = T/2. In this case T = 2L and the step scaling technique (see eq.(4)) 
connects xq = L m i n , where the meson mass has been defined on the smallest 
volume, with x$ = L max , where one expects to be free from both sources of 
finite volume effects. 

The finite volume effects due to the excited states can be easily understood 
by means of a two mass model for the correlation 



Here A(L) and B(L) are coefficients, M(L) is the meson mass and 8M(L) 
is the mass shift. All these quantities depend upon the physical extension of 




(16) 



fx(x ,L) = A(L) e - M ^ X0 + B(L) e ~WL)+SA4(L)] Xo 



(17) 
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the volume L but, in order to isolate the effects due to the excited states, we 
consider M and 8M constants. In this simplified model the effective mass of 
eq. (16) takes the form 

r 1 + B(L) e -8M(x -a) 

aM x (x ) =aM+- In I ^ \ (18) 

As in numerical simulations, we take xq = L. The functional dependence upon 
the volume of the ratio B(L)/A(L) can be inferred from the data on the meson 
decay constants given in [8,9]. A parametrization that fits well the data is given 
by 

A(L) = a (l + ^ B(L) = -0.5 (19) 

In Fig. [1] are shown the plots of the effective mass step scaling function, as 
derived from this model, defined by 

as functions of the volume and at different values of the meson masses. At the 
volume L = 0.4 fm the step scaling functions are smaller than one while, the 
pattern is reversed at L = 0.8 fm. Independently from the volume, the heavier 
mesons are closer to unity than the lighter ones. 

The qualitative predictions of this simple model well reproduces the behavior 
of the numerically measured step scaling functions, shown in Fig. [3] for L = 
0.4 fm and in Fig. [6] for L = 0.8 fm. 



4 Step scaling functions and HQET 

In this section we enter into the details of the application of the SSM, 
illustrated in sec. 2, to the particular case of the heavy-light meson masses 
computation. We introduce the step scaling functions 



rr (T r n RGI ™ RGI } 

(7p \L/, Tfli , 77l 2 J 



M P 


[m RGI , 


m RGI ^j 


2L 


M P 


(m RGI 


m RGI ) 


L 


My 


(mf GI , 


m RGI ^j 


2L 


My 


(m RGI 


m RGI ) 


L 



(t RGI RGl\ _ V 1 ' 2 J 2L /oi\ 



for the pseudoscalar and vector meson masses respectively. 

To validate the expansion of eq. (3), we can make use of the HQET predictions 

on the heavy-light meson masses. In the infinite volume the pseudoscalar and 
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vector meson masses have the following expansion in terms of the heavy-quark 
mass: 

M x (m h , mi ) = m h + K(m l ) + ax{m <^ + ,.. (22) 



where X e {P, V}. Assuming the contribution of the l/m\ corrections to be 
negligible, at finite volume one has 

M x (m h , m h L)=m h + k x (m h L) + ax ^ L ^ + . . . (23) 



where A x (mi, L) depends upon the spin of the meson state because of the 
contamination of the excited states to the finite volume correlations. Using 
eqs. (21) and (23) we obtain the HQET predictions for the step scaling func- 
tions of the heavy-light meson masses 

E x (L,m h , mi ) = 1+ xK - } + * V ; +... 24 

m h mi 



This result requires some considerations. First we want to stress that in the 
infinite heavy-quark mass limit the step scaling functions have to be exactly 
equal to one, T<x(L,mi,mh — > 00) = 1 1 . This represents a strong constraint 
for the fits of the heavy-quark mass dependence of the step scaling functions. 

The second observation concerns the number of terms to be considered in 
eq. (24). At order 0(1 /nth) one has 

Y,%\m h L) = K x (m h 2L) - A x (m h L) (25) 



corresponding to the static approximation in eq. (23). By increasing the phys- 
ical volume L, the difference between A x (mi,2L) and A x (mi,L) decreases 
because the two quantity have to be equal in the infinite volume limit, making 
the heavy-quark mass expansion of the finite volume effects rapidly conver- 
gent. 

The same arguments apply to the coefficient E^(mj, L) that has to be consid- 
ered when in the expansion of the meson masses, eq. (23), the order 0(1 /mh) 
is taken into account. 

In our calculation we will perform the fits of the step scaling functions consid- 
ering the 0(1 /ml) term, beyond the so called static evolution (SE) that 
retains only. We will also report for comparison the fits for the SE that 
are anyway compatible. 



We thank A. Kronfeld and R. Sommer for having pointed out this property of 
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a 


L/a 




k 


m RGi (GeV) 








0.120081 


7.14(8) 








0.120yoo 


6.63(7) 








1 26050 


4.024(44) 


6.963 


16 


0.134827(6) 


0.131082 


1.696(19) 








0.131314 


1.591(18) 








0.134526 


0.1381(30) 








0.134614 


0.0978(28) 








0.134702 


0.0574(28) 








0. 1Z41 ( D 


7.11(8) 








f\ "1 O A O A A 

0.124844 


6.61(20) 








0.128440 


4 ni8('44 , | 


7.300 


24 


0.134235(3) 


0.131800 


1.695(19) 








0.131950 


1.592(18) 








0.134041 


0.1374(27) 








0.134098 


0.0971(24) 








0.134155 


0.0567(24) 








0.126352 


7.10(8) 








0.126866 


6.60(7) 








0.129585 


4.016(44) 


7.548 


32 


0.133838(2) 


0.132053 


1.698(19) 








0.132162 


1.595(18) 








0.133690 


0.1422(27) 








0.133732 


0.1021(25) 








0.133773 


0.0618(23) 



Table 1 

Simulation parameters at Lq = 0.4 fm. The RGI quark masses are obtained using 
eq. (12). 

Similar arguments could be repeated for the heavy-heavy step scaling 
functions using NRQCD. In this case one cannot predict the value of the 
step scaling functions at rrih — > oo. With the three simulated values of the 
heavy quark masses we extrapolate our numerical data with a simplified NLO 
formula, i.e. a linear dependence upon the heavy-quark mass, given the little 
sensitivity to possible slowly varying logarithmic terms. 



5 Numerical results 

This section contains the simulation parameters and the numerical results 
for each step of the calculation. 
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Fig. 2. Continuum extrapolation on the small volume, Lq, of the mass Mp(Lq) 
of the pseudoscalar heavy-heavy meson corresponding to the heavy quark of mass 
m = 7.10 GeV. The different values of the meson mass correspond to different 
definitions of the RGI quark masses given in equations (12), (13) and (14). Units are 
in GeV. Similar plots can be obtained, from the data reported in Table [6,7,8], for 
the other combinations of quark masses used in our simulations, also in the case of 

My(L ). 

5.1 The small volume: Lq = 0.4 fm 



The physical extension of the small volume has been chosen in order to 
properly account the dynamics of quarks with masses in the region of the the 
physical 5-quark; we have fixed it to be L = 0.4 fm. In order to have a contin- 
uum extrapolation of the numerical results, this volume has been simulated 
using three different discretization, 32 x 16 3 , 48 x 24 3 and 64 x 32 3 . Using 
the results reported in [25] we have fixed two heavy quark masses with values 
interpolating the 6-quark mass. In the renormalization group invariant (RGI) 
scheme fixed by the conventions reported in [22,16], these masses are 7.10 
GeV and 6.60 GeV respectively. Two other heavy quarks have been simulated 
in order to interpolate the mass region where, using the results reported in 
[26,25], one expects to find the physical c-quark, i.e. 1.70 GeV and 1.60 GeV 
respectively. The charm quark mass is affordable also by a direct computation: 
here the calculation via our step scaling method represents a check of the pro- 
cedure. An additional heavy quark has been simulated with mass 4.00 GeV. 
Three light quark have been simulated with masses of 0.14 GeV, 0.10 GeV and 
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Lo/a k c 


k 


m RGi (GeV) 






0.120674 


3.543(39) 






0.122220 


3.114(34) 






0.126937 


1.927(21) 


6.420 


8 0.135703(9) 


0.134304 


0.3007(36) 






0.134770 


0.2003(28) 






0.135221 


0.1028(21) 






0.1249 


3.542(39) 






0.1260 


3.136(34) 






0.1293 


1.979(22) 


6.737 


12 0.135235(5) 


0.1343 


0.3127(38) 






0.1346 


0.2090(28) 






0.1349 


0.1080(21) 






0.127074 


3.549(39) 






0.127913 


3.153(35) 






0.130409 


2.003(22) 


6.963 


16 0.134832(4) 


0.134145 


0.3134(38) 






0.134369 


0.2112(28) 






0.134593 


0.1086(20) 



Table 2 

Simulation parameters for the first evolution step Lq — > L\ = 0.8 fm. The RGI quark 
masses are obtained using eq. (12). 

0.06 GeV. Using the accurate determination of the RGI strange quark mass 
given in [27] we have fixed one of the simulated light quarks to be the physi- 
cal s. All the parameters of the three different simulations are summarized in 
Table [1]. The numerical results of the pseudoscalar and vector meson masses, 
M P and My, are shown in Table [6,7,8] both for the heavy-heavy and for the 
heavy-light quark anti-quark pairs. 

In order to obtain physical predictions from the simulated data on this 
finite volume we need to extrapolate our numerical results to the continuum. 
As already mentioned, we have obtained different set of data by using the 
different definitions of the RGI quark masses given in the equations (12), (13) 
and (14). The continuum results are thus obtained trough a combined fit of all 
the set of data, linear in (a/r ) 2 , as shown in Fig. [2] in the case of the mass 
of the pseudoscalar heavy-heavy meson corresponding to the heavy quark of 
mass m RGI = 7.10 GeV. We obtain a global x 2 /dof = 5.40 to be compared 
with the x 2 s of each individual definition listed in the figure. We also show 
the points at the largest lattice spacing not included in the fit. 
The errors included in the evaluation of the y 2 are statistical only. These 
are calculated by a jackknife procedure, also in the case of the step scaling 
functions. The systematics due to the uncertainty on the lattice spacing has 
been estimated by repeating the fit using the different values of the scale 
allowed by the uncertainties quoted in [10,11,12] and considering the spread 
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of the results. The same procedure has been used for the systematics due to 
the uncertainties on the renormalization constants. The resulting 2% percent 
for the renormalization constants and 1% percent for the scale are summed in 
quadrature and added to the statistical errors. 

5.2 The first evolution step: L i— > L\ = 0.8 fm 

The finite volume effects on the quantities calculated at Lq, are measured 
doubling the volume, L\ = 0.8 fm, by using the step scaling functions of 
eqs. (21). 

In order to have results in the continuum, also in the case of the step scaling 
functions three different discretizations of Lq have been used, i.e 16 x 8 3 , 
24 x 12 3 and 32 x 16 3 . The volume L\ has been simulated starting from the 
discretizations of L , fixing the value of the bare coupling and doubling the 
number of lattice points in each direction. 
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Fig. 3. The figure shows the pseudoscalar step scaling functions £p as functions of 
rrii GI , for the simulation of the first evolution step corresponding to (3 = 6.963. The 
different sets of data correspond to the values of . As can be seen the step scaling 
functions approach a plateau for high values of mf G1 . Similar plots can be obtained 
using the data of Table [9,10,11] for the other values of the bare coupling and for the 
vector mesons step scaling functions £y. 

The simulated quark masses have been halved with respect to the masses 
simulated on the small volume in order to have the same order of discretization 
effects proportional to am. The set of parameters for the simulations of this 
evolution step is reported in Table [2] and the numerical results are given in 
Table [9,10,11]. 

The step scaling functions of the pseudoscalar mesons at (3 — 6.963 are 
plotted, at fixed mf G/ , as functions of mf G/ in Fig. [3]. The value of the step 
scaling functions for the s quark are obtained trough linear interpolation. 

Both Sp and Sy are almost flat in a region of heavy quark masses start- 
ing around the charm mass. The hypothesis of low sensitivity upon the high- 
energy scale is thus verified. 

In Fig. [4] are reported the results of the continuum extrapolation of the step 
scaling function, Sp(L ), of the pseudoscalar meson corresponding the heavy 
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Fig. 4. Continuum extrapolation on the first evolution step, Lq h-> Lj, of the step 
scaling function, Ep(Ln), of the pseudoscalar meson corresponding to the heavy quark 
of mass mf G/ = 3.55 GeV. The two sets of data are obtained using the two definitions 
of RGI quark masses of equations (12) and (14). Units are in GeV. Similar plots can 
be obtained, from the data reported in Table [9,10,11], for the other combinations of 
quark masses used in our simulations, also in the case of Sy(Lo). 
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Fig. 5. The figure shows the continuum extrapolated step scaling functions op(Lo) as 
functions of l/mi G1 . The heavy extrapolations are shown only for the heavy-strange 
(Hs) set of data. Similar plots can be obtained using the data of Table [9,10,11] for 
the vector mesons step scaling functions ay(Lo). 

quark of mass mf G/ = 3.55 GeV. 

The residual heavy mass dependence of the continuum extrapolated step scal- 
ing functions is very mild, as shown in Fig. [5] in the plot of op as a function 
of the inverse quark mass. 

In the heavy-light already discussed in sec. 4, we extract the values of 

the step scaling functions at the values of the heavy-quark masses simulated 
on the small volume by interpolation between the numerical results and the 
theoretical point at m^ GI — > oo using both a linear fit (SE) and a quadratic 
fit (RE). In the heavy-heavy case the results are linearly extrapolated. 



5.3 The second evolution step: L\ i— > L2 — 1.6 fm 



In order to have the results on a physical volume, L 2 = 1.6 fm, a second 
evolution step it is required. This is done computing the meson mass step 
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Fig. 6. The figure shows the pseudoscalar step scaling functions £p as functions of 
mf G/ , for the simulation of the second evolution step corresponding to f3 = 6.420. The 
different sets of data correspond to the values of . Similar graphs can be obtained 
using the data of Table [12,13,14] for the other values of the bare coupling and for 
the vector mesons step scaling functions £y. 
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Fig. 7. Continuum extrapolation on the second evolution step, L\ i— » L2, of the step 
scaling function, Sp(Li), of the pseudoscalar meson corresponding to the heavy quark 
of mass — 2.00 GeV. The two sets of data are obtained using the two definitions 

of RGI quark masses of equations (12) and (14). Units are in GeV. Similar plots can be 
obtained, from the data reported in Table [12,13,14], for the other combinations of 
quark masses used in our simulations, also in the case of Xy(Li). 
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Fig. 8. The figure shows the continuum extrapolated step scaling functions ap(Li) as 
functions of l/m^ 1 . The heavy extrapolations are shown only for the heavy-strange 
(Hs) set of data. Similar plots can be obtained using the data of Table [12,13,14] for 
the vector mesons step scaling functions oy(Li). 
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m RGI ( GeV) 






0.118128 


2.012(22) 






0.121012 


1.551(17) 






0.122513 


1.337(15) 


5.960 


8 0.13490(4) 


0.131457 


0.3154(36) 






0.132335 


0.2322(28) 






0.133226 


0.1466(44) 






0.124090 


1.984(22) 






0.126198 


1.584(17) 






0.127280 


1.389(15) 


6.211 


12 0.135831(8) 


0.133574 


0.3493(39) 






0.134177 


0.2550(29) 






0.134786 


0.1510(19) 






0.126996 


1.933(21) 






0.128646 


1.547(17) 






0.129487 


1.355(14) 


6.420 


16 0.135734(5) 


0.134318 


0.3016(34) 






0.134775 


0.2038(24) 






0.135235 


0.1055(15) 



Table 3 

Simulation parameters for the first evolution step L\ — > Li = 1.6 fm. The RGI quark 
masses are obtained using eq. (12). 



scaling functions of eq. (21) at L — Li, by the procedure outlined in the 
previous section. The parameters of the simulations are given in Table [3] and 
the results are in Table [12,13,14]. 

Also in this case, the values of the simulated quark masses have been 
halved with respect to the previous step, owing to the lower values of the 
simulation cutoffs. In Fig. [6] we show the pseudoscalar step scaling functions at 
(3 = 6.420 and in Fig. [8] the residual heavy-quark mass dependence with the SE 
and RE fits. 

Fig. [7] shows the continuum extrapolation of the step scaling function, Sp(Li), 
of the pseudoscalar meson corresponding to the heavy quark of mass mf 01 = 2.00 
GeV. 

The contributions of the excited states, predicted from the simple model of 
eq. (20) and present in our numerical data, should disappear in a third evolution 
step. A check supporting this hypothesis is that on the larger volumes used in the 
simulation of this evolution step, i.e. L2 at (3 = 6.420 and (3 = 6.211, the values of 
the meson masses defined at xq = T/2 coincide, within the errors, with the values 
coming from a single mass fit to the correlations. 
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hh 
hs 
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6.18(14) 
3.55(8) 
3.46(10) 


6.20(14) 
3.62(8) 
3.53(11) 


1.70 


hh 
hs 
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3.15(7) 
1.97(5) 
1.88(6) 


3.21(7) 
2.09(5) 
2.00(6) 


1.60 


hh 
hs 
hu 


3.01(7) 
1.90(5) 
1.80(5) 


3.07(7) 
2.01(5) 
1.93(6) 



Table 4 

Meson masses in the infinite volume limit. Units are in GeV. 




Fig. 9. Chiral extrapolation of the continuum results for the pseudoscalar meson cor- 
responding to m^ GI = 7.1. Units are in GeV. 

6 Physical results 



In this section we combine the results of the small volume with the results of 
the step scaling functions to obtain, according to eq. (4), the physical numbers. 
In Table [4] we give the pseudoscalar and vector meson masses corresponding to 
the heavy quarks simulated on the small volume and, in the heavy-light case, to 
strange and up light quarks. The results corresponding to the up quark have been 
obtained by extrapolating the continuum data in the large volume from masses 
around the strange region. In Figure [9] we show the extrapolation for the pseu- 
doscalar meson corresponding to m? =7.1 GeV. Having simulated three light 
quark on the largest volume we have fitted the results linearly without trying com- 
plicated functional forms requiring more than two parameters. Comparing these 
results with the experimental determinations of the same quantities [25] we obtain 
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different determinations of the 6-quark mass, depending upon the physical state 
used as experimental input. The results are summarized in Table [5]. 







m h 






from P 




from V 


lb 


6.44(14) 




6.57(14) 


is 


6.91(16) 




6.93(16) 


bu 


6.90(20) 




6.91(20) 


cc 


1.603(35) 




1.642(36) 


cs 


1.692(38) 




1.741(39) 


cu 


1.690(50) 




1.712(51) 



Table 5 

Determinations of the heavy quark masses from the heavy-heavy and from the heavy- 
light states. The errors include our estimate of the systematics. Units are in GeV. 

Within the quenched approximation, the determinations of the quark masses 
coming from the heavy-heavy or from the heavy-light spectrum in principle differ 
because the theory does not account for the fermion loops. We obtain two deter- 
minations that are marginally compatible within the errors and that might suggest 
the need for a tiny unquenching effect. The good agreement between the determi- 
nations of the quark masses coming from the heavy-up and heavy-strange sets of 
data make us confident on our chiral extrapolations. Of course the heavy-strange 
and the heavy-heavy case are not extrapolated at all and we therefore use these 
results only to determine the heavy-quark masses. The numbers we quote as final 
results are obtained by averaging the four results in the first two rows of the two 
sets of Table [5] for the heavy-heavy and the heavy-strange cases and by keeping 
the typical error of a single case: 

m^ GI = 6.73(16) GeV mf^(mf^) = 4.33(10) GeV (26) 

for the 6-quark and 

m f GI = 1.681(36) GeV mf^mf^) = 1.319(28) GeV (27) 

for the charm. The latter results compare favorably with the results of the direct 
computations [26,25]. 

As already explained in section 5.1, our error estimate includes both the statis- 
tical error from the Monte Carlo simulation as well as the systematic error coming 
from the uncertainty on the lattice spacing corresponding at a given (3 value and 
to the uncertainty on the renormalization constants in eqs (10) and (14). The final 
errors on the continuum quantities, of the order of 2% percent for the renormaliza- 
tion constants and of about 1% percent for the scale, are summed in quadrature 
and added to the statistical errors. The evolution to the MS scheme has been done 
using four-loop renormalization group equations [28,29,30]. 

Using these determinations of the quark masses we can calculate the mass of the 
B c meson from our set of correlations suitably interpolated. Following this procedure 
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we obtain a determination that is in very good agreement with the experimental 
measurement of the meson mass [25] : 

M Bc = 6.46(15) GeV M^f = 6.40(39) (13) GeV (28) 

The same analysis has been done also in the case of the decay constants and in [9] we 
give the first determination of fs c from lattice QCD in the quenched approximation. 

7 Conclusions 

We have shown the effectiveness of the step scaling method in the computation 
of the heavy-light and heavy-heavy meson spectrum. We have performed the cal- 
culations for the pseudoscalar and vector mesons in the continuum limit of lattice 
QCD in the quenched approximation and, by comparison with the experimental de- 
terminations, we have extracted the masses of the bottom and of the charm quarks. 
Our results represent the first determination of the bottom quark mass in the con- 
tinuum limit of quenched lattice QCD. 

Using the results on the quark masses we give an independent theoretical cal- 
culation of the B c meson mass that matches very well the experimental result. 
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Table 6 

Values of the pseudoscalar, Mp, and vector, My 
simulations on the smallest volume Lq = 0.4 fm at 



2.771(9) 
1.888(12) 
1.864(12) 
1.840(12) 



meson masses resulting from the 
P = 6.963. Units are in GeV. 
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Table 7 

Values of the pseudoscalar, Mp, and vector, My 
simulations on the smallest volume Lq = 0.4 fm at 



2.813(13) 
1.934(17) 
1.911(17) 
1.888(18) 



meson masses resulting from the 
P = 7.300. Units are in GeV. 
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Table 8 

Values of the pseudoscalar, Mp, and vector, My 
simulations on the smallest volume Lq = 0.4 fm at 



2.802(13) 
1.924(17) 
1.901(18) 
1.879(18) 



meson masses resulting from the 
P = 7.548. Units are in GeV. 
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Table 9 

Values of the step scaling functions for the pseudoscalar, £p, and vector, Sy, mi 
masses resulting from the simulation at (3 = 6.420 of the first evolution step Lq — > 
Units are in GeV. 
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0.1080(21) 



0.1080(21) 
0.2090(28) 
0.3127(38) 
1.979(22) 
3.136(34) 
3.542(39) 



0.907(25) 
0.935(21) 
0.951(19) 
1.009(7) 
1.014(5) 
1.014(5) 



0.789(19) 
0.810(17) 
0.830(16) 
0.976(7) 
0.995(5) 
0.998(5) 



0.2090(28) 0.2090(28) 0.950(18) 0.828(16) 

0.3127(38) 0.960(16) 0.846(14) 

1.979(22) 1.011(6) 0.980(7) 

3.136(34) 1.0155(47) 0.998(5) 

3.542(39) 1.0158(43) 1.0009(46) 



6.737 


0.3127(38) 


0.3127(38) 
1.979(22) 
3.136(34) 
3.542(39) 


0.967(14) 
1.014(6) 
1.0176(44) 
1.0178(40) 


0.862(13) 
0.985(6) 
1.0013(47) 
1.0039(43) 




1.979(22) 


1.979(22) 
3.136(34) 
3.542(39) 


1.0376(30) 
1.0375(23) 
1.0367(21) 


1.0279(33) 
1.0317(25) 
1.0317(23) 




3.136(34) 


3.136(34) 
3.542(39) 


1.0372(18) 
1.0365(16) 


1.0338(19) 
1.0336(18) 




3.542(39) 


3.542(39) 


1.0358(15) 


1.0333(17) 



Table 10 

Values of the step scaling functions for the pseudoscalar, £p, and vector, £y, mi 
masses resulting from the simulation at (3 = 6.737 of the first evolution step Lq — > 
Units are in GeV. 
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0.1086(20) 



0.1086(20) 
0.2112(28) 
0.3134(38) 
2.003(22) 
3.153(35) 
3.549(39) 



0.914(28) 
0.939(24) 
0.953(21) 
1.012(8) 
1.016(6) 
1.016(5) 



0.783(24) 
0.807(21) 
0.830(19) 
0.981(8) 
0.999(6) 
1.002(6) 



0.2112(28) 
0.3134(38) 
2.003(22) 
3.153(35) 
3.549(39) 



0.951(21) 
0.960(18) 
1.013(7) 
1.017(5) 
1.017(5) 



0.828(19) 
0.847(17) 
0.984(8) 
1.001(6) 
1.003(5) 



6.963 


0.3134(38) 


0.3134(38) 
2.003(22) 
3.153(35) 
3.549(39) 


0.967(16) 
1.015(7) 
1.018(5) 

1.0181(46) 


0.864(16) 
0.989(7) 
1.004(5) 
1.006(5) 




2.003(22) 


2.003(22) 
3.153(35) 
3.549(39) 


1.0365(35) 
1.0357(27) 
1.0348(25) 


1.0281(39) 
1.0309(30) 
1.0307(28) 




3.153(35) 


3.153(35) 
3.549(39) 


1.0349(21) 
1.0341(20) 


1.0321(23) 
1.0317(22) 




3.549(39) 


3.549(39) 


1.0333(19) 


1.0313(20) 



Table 11 

Values of the step scaling functions for the pseudoscalar, £p, and vector, £y, mi 
masses resulting from the simulation at (3 = 6.963 of the first evolution step Lq — > 
Units are in GeV. 
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0.1466(44) 



0.1466(44) 
0.2322(28) 
0.3154(36) 
1.337(15) 
1.551(17) 
2.012(22) 



1.269(13) 
1.244(11) 
1.226(9) 
1.1406(40) 
1.1332(37) 
1.1214(32) 



1.308(15) 
1.286(12) 
1.267(10) 
1.1636(45) 
1.1539(42) 
1.1383(36) 



0.2322(28) 
0.3154(36) 
1.337(15) 
1.551(17) 
2.012(22) 



1.224(9) 
1.209(8) 
1.1337(35) 
1.1268(32) 
1.1158(29) 



1.266(10) 
1.249(9) 
1.1557(40) 
1.1467(37) 
1.1322(32) 



5.960 


0.3154(36) 


0.3154(36) 
1.337(15) 
1.551(17) 
2.012(22) 


1.196(7) 
1.1278(32) 
1.1213(30) 
1.1110(27) 


1.235(8) 
1.1489(37) 
1.1405(34) 
1.1269(30) 




1.337(15) 


1.337(15) 
1.551(17) 
2.012(22) 


1.0894(19) 
1.0853(18) 
1.0785(16) 


1.029(21) 
1.0978(20) 
1.0895(17) 




1.551(17) 


1.551(17) 
2.012(22) 


1.0814(17) 
1.0749(15) 


1.0931(18) 
1.0852(16) 




2.012(22) 


2.012(22) 


1.0689(13) 


1.0782(14) 



Table 12 

Values of the step scaling functions for the pseudoscalar, Sp, and vector, Sy, mi 
masses resulting from the simulation at (3 = 5.960 of the second evolution step L\ — > 
Units are in GeV. 
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0.1510(19) 



m% GI 

0.1510(19) 
0.2550(29) 
0.3493(39) 
1.389(15) 
1.584(17) 
1.984(22) 



1.156(14) 
1.151(11) 
1.146(10) 
1.1035(46) 
1.1098(42) 
1.0899(37) 



1.255(18) 
1.235(15) 
1.219(12) 
1.131(6) 
1.123(5) 
1.1095(44) 



0.2550(29) 0.2550(29) 1.147(10) 1.218(12) 

0.3493(39) 1.142(8) 1.204(10) 

1.389(15) 1.1010(40) 1.1258(48) 

1.584(17) 1.0961(37) 1.1181(44) 

1.984(22) 1.0879(33) 1.1056(39) 



6.211 


0.3493(39) 


0.3493(39) 
1.389(15) 
1.584(17) 
1.984(22) 


1.138(7) 
1.0983(37) 
1.0936(34) 
1.0857(30) 


1.193(9) 
1.1210(44) 
1.1137(40) 
1.1020(35) 




1.389(15) 


1.389(15) 
1.584(17) 
1.984(22) 


1.0727(20) 
1.0695(19) 
1.0640(17) 


1.0842(24) 
1.0799(22) 
1.0728(20) 




1.584(17) 


1.584(17) 
1.984(22) 


1.0664(18) 
1.0611(16) 


1.0759(21) 
1.0692(18) 




1.984(22) 


1.984(22) 


1.0564(14) 


1.0633(16) 



Table 13 

Values of the step scaling functions for the pseudoscalar, Sp, and vector, Sy, mi 
masses resulting from the simulation at (3 = 6.211 of the second evolution step L\ — > 
Units are in GeV. 
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0.1055(15) 
0.2038(24) 
0.3016(34) 
1.355(14) 
1.547(17) 
1.933(21) 



1.167(27) 
1.162(22) 
1.156(18) 
1.104(7) 
1.098(7) 
1.088(6) 



1.259(35) 
1.235(27) 
1.218(22) 
1.125(9) 
1.117(8) 
1.103(7) 



0.2038(24) 0.2038(24) 1.155(18) 1.217(21) 

0.3016(34) 1.150(15) 1.203(18) 

1.355(14) 1.101(6) 1.121(8) 

1.547(17) 1.095(6) 1.113(7) 

1.933(21) 1.086(5) 1.000(6) 



6.420 


0.3016(34) 


0.3016(34) 
1.355(14) 
1.547(17) 
1.933(21) 


1.144(13) 
1.098(6) 
1.093(5) 

1.0839(46) 


1.192(16) 
1.118(7) 
1.110(6) 
1.098(5) 




1.355(14) 


1.355(14) 
1.547(17) 
1.933(21) 


1.0711(31) 
1.0676(29) 
1.0617(25) 


1.0821(37) 
1.0777(34) 
1.0702(30) 




1.547(17) 


1.547(17) 
1.933(21) 


1.0643(27) 
1.0587(24) 


1.0735(31) 
1.0665(27) 




1.933(21) 


1.933(21) 


1.0538(21) 


1.0604(24) 



Table 14 

Values of the step scaling functions for the pseudoscalar, Sp, and vector, Sy, mi 
masses resulting from the simulation at (3 = 6.420 of the second evolution step L\ — > 
Units are in GeV. 
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